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I  .   INTRODUCTION 

One  of  the  dominant  features  of  ocean  waves  is  the  occur- 
rence of  wave  groups.   The  existence  and  properties  of  wave 
groups  have  become  increasingly  important  for  economic 
reasons.   Groups  of  large  waves  have  been  damaging  to  coastal 
structures,  both  man-made  and  natural.   Surf  zone  water  level 
fluctuations  (surf  beat)  and  current  variations  are  also  re- 
lated to  wave  groups.   Induced  seiching  in  harbors  may  occur, 
which  may  lead  to  stress  on  moored  craft,  and  in  extreme 
cases  to  flooding  of  adjacent  land  areas.   Wave  groups  in- 
fluence the  response  of  ships  at  sea.   The  literature  devoted 
to  this  topic  has  been  fairly  sparse  in  the  past,  but  an  in- 
creasing interest  is  becoming  apparent.   Papers  dealing 
specifically  with  the  definition,  occurrence,  and  properties 
of  wave  groups  include  Goda  (1970)  ,  Thompson  (1972)  ,  Ewing 
(1973),  Smith  (1974),  Siefert  (1976),  and  Chou  (1978). 

Properties  that  are  considered  to  define  "wave  groups" 
are:  (1)  groups  consist  of  consecutive  large  waves,  and  (2) 
the  waves  contained  within  a  group  are  nearly  periodic.   In 
this  study  a  simple  method  is  developed,  which  can  be  per- 
formed with  computer  assistance,  for  the  identification  of 
wave  groups  in  a  wave  record.   The  method  involves  analysis 
of  wave  records  with  a  type  of  rms  smoothing  filter.   The 
smoothed  record  is  then  used  to  identify  groups  and  determine 
appropriate  measures  of  wave  group  properties.   The  smoothing 
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mechanism  is  a  sliding  short-term  variance.   Relationships 
are  then  probed  between  wave  group  measures  and  wave  record 
measures,  and  lastly,  the  interrelationships  of  the  differ- 
ent wave  group  measures  are  examined.   The  analyzed  wave  rec- 
ords contain  ocean  swell  representing  a  wide  range  of  wave 
steepness,  and  were  recbrded  on  the  Southern  California 
coast . 
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II.   WAVE  GROUP  DETERMINATION 

A.   FOUNDATIONAL  CONSIDERATIONS 

There  are  two  fundamental  approaches  to  the  analysis  of 
wave  records,  an  analog  method  which  is  concerned  with  dis- 
tributions of  wave  maxima  and  minima,  i.e.,  crests  and 
troughs  that  occur  at  irregular  time  intervals,  and  a  digital 
method  which  allows  examination  of  wave  records  at  equally 
spaced  time  intervals  without  regard  to  crests  and  troughs. 
The  latter  approach  is  used  in  this  investigation  to  define 
and  identify  wave  groups,  and  to  measure  their  energy  content. 

To  process  wave  data  in  a  search  for  wave  group  proper- 
ties through  the  use  of  a  computer  requires  that  some  wave 
record  characteristics  be  known.   For  ease  of  computation, 
and  in  consonance  with  past  studies,  the  mean  water  level 
will  be  considered  to  be  the  average  of  the  water-level  devia- 
tion from  a  standard  baseline  (commonly  the  chart  bottom  in 
an  analog  record)  for  an  entire  record  (in  this  study  approxi- 
mately 17-18  minutes).   The  waves  are  assumed  to  approximate 
a  symmetrical  distribution  about  this  level.   For  the  purposes 
of  this  study,  an  individual  wave  is  defined  as  beginning  at 
one  upcrossing  of  the  mean  water  level  by  the  surface  of  the 
water,  and  ending  at  the  next  adjacent  upcrossing  in  time. 

The  procedure  for  the  identification  of  a  wave  group  is 
developed  from  the  concepts  and  interrelationships  of  the 
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variance,  energy,  and  height  parameters  of  a  wave  record. 
The  sample  variance,  V,  of  a  finite  length  wave  record  may  be 
determined  from  the  following  equation: 


V  *  nTT  £   <xi  "  x> 
i  =  l 


-.  2 


(1) 


where  n  is  the  number  of  equally  spaced  samples  in  time,  x. 
is  the  water  surface  elevation  at  time  i,  and  x  is  the  sample 
mean  water  level  found  by 


n  *-rf    1 


i=l 

This  study  uses  the  concept  that  the  energy  represented 
in  a  wave  record,  the  area  under  a  spectral  density  curve,  is 
proportional  to  the  variance  of  that  record.   The  energy 

(average  energy  per  unit  area  of  sea  surface)  of  a  simple 

2 
sinusoidal  wave  can  be  written  as  E  =  pgH  /8,  where  H  is  the 

wave  height  from  crest  to  trough.   If  the  sea  surface  is 
thought  of  as  the  sum  of  an  infinite  number  of  small  sinus- 
oidal waves,  then  the  total  wave  energy  present  may  be  con- 
sidered to  oe  the  sum  of  the  energies  contained  in  the  indi- 
vidual waveforms,  which  is  proportional  to  the  sum  of  the 
squares  of  the  heights  of  those  waveforms  (Michel,  1968)  . 
The  variance  is  then  related  to  the  sum  of  the  squares  of  the 
wave  heights . 
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This  study  assumes  a  Rayleigh  distribution  of  ocean  wave 
heights  (Longuet-Higgins ,  1952).   Its  probability  density 
function  can  be  written  in  terms  of  the  variance,  V, 

p(H)  =|e"H2/2V  (2) 

This  distribution  describes  the  frequency  of  occurrence  of 
wave  height  H  occurring  in  a  finite  set  of  recorded  waves. 
Thus,  through  the  use  of  the  Rayleigh  formula,  if  the  vari- 
ance or  energy  of  the  record  is  known,  the  distribution  of 
wave  heights  may  be  found.   Certain  height  parameters  may 
therefore  be  found  directly  from  the  energy,  or  variance,  of 
a  record  using  equations  of  the  form  H  =  k/v,  where  k  is  a 
constant  (Michel,  1968). 

The  significant  height,  H  ,  the  average  height  of  the 
highest  one-third  waves,  and  the  corresponding  variance  will 
be  used  herein  for  the  development  of  a  method  of  defining 
wave  groups.   These  parameters  are  related  by  H   =  4/v. 

B.   INDIVIDUAL  WAVE  GROUP  DETERMINATION 
1 .   Definition  of  a  Wave  Group 

Wave  groups  are  noticeable  to  a  casual  observer  be- 
cause the  heights  of  the  successive  waves  comprising  the 
group  differ  from  a  succession  of  random  wave  heights.   The 
short-term  variance  of  a  digitized  wave  record  over  an  inter- 
val, such  as  the  time  a  wave  group  is  present,  is  thus  larger 
than  the  variance  of  the  wave  record  as  a  whole. 

With  this  concept  in  mind,  individual  wave  records 
could  be  examined  by  looking  through  a  window  of  limited 
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duration  r  to  obtain  a  short  sample  of  the  record,  as  in  Fig- 
ure 1.   The  variance  of  this  sample,  V.,  may  then  be  deter- 
mined . for  t.  =  iAt  by  the  use  of 
1 


1  + 


V, 


n-1   2-f 


3  =  1- 


(x,  - 


x)2 


(3) 


where  t.  is  the  midpoint  of  the  window,  At  is  the  sampling  in- 
terval of  the  record  (At  equals  one  second  for  this  study) , 
and  x  is  the  mean  of  the  entire  record.   The  number  of  samples 
(n)  in  the  window  and  its  width  (r)  are  related  by  the  equa- 
tion r  =  nAt.   In  Figure  1,  a  =  (i  -  — )  At  and  b  =  (i  +  — )  At. 
The  variance  is  converted  into  a  significant  height.   This 

height,  H   ,.  is  plotted  on  a  graph  according  to  the  midpoint 
S  G 

of  the  window,  (a+b)/2  (Figure  1).   The  significant  height  of 

the  entire  observed  record,  H   ,  is  also  plotted.   A  height 

S  R 

parameter  (proportional  to  the  square  root  of  the  variance) 
was  used  in  this  study  rather  than  the  variance  for  convenience 
in  graphical  scaling  and  the  significant  height  was  chosen 
as  this  parameter  because  of  familiarity. 

In  smoothing  a  record  for  wave  group  characteristics 
the  window  is  moved  through  the  entire  record  in  time  incre- 
ments corresponding  to  the  sampling  rate  of  the  record,  and 
the  significant  height  values  plotted.   Interconnection  of 

these  values,  H   ,  produces  a  waveform,  called  the  group  wave- 
S  G 

form,  which  may  then  be  compared  to  the  value  of  the  signifi- 
cant height  of  the  entire  record,  H    (Figure  2).   A  portion 

S  R. 
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of  this  new  record  where  the  H    waveform  is  greater  than  H 

S  G  jR 

is  a  region  where  the  variance  over  a  short  time  span  is 
greater  than  the  variance  of  the  record  as  a  whole.   This  re- 
gion of  higher  values  is  considered  to  identify  a  wave  group. 

The  term  "wave  group"  implies  a  series  of  whole  waves, 
and  it  can  be  seen  that  this  region  (where  H    >  H   )  does 

o  G       dK 

not  necessarily  include  an  integral  number  of  waves,  nor  do 
its  boundaries  coincide  with  some  distinctive  part  of  the  ob- 
served waveform  such  as  the  crest,  trough,  zero-upcrossing , 
or  zero-downcross ing  point.   The  determination  of  the  wave 
group  boundaries  adopted  in  this  study  is  illustrated  in  Fig- 
ure 3.   The  group  waveform  computed  using  a  moving  time  window 

is  seen  to  intersect  the  value  of  H  „  at  times  t   and  t,  ,  dur- 

SR  a       b 

ing  which  H    >  H   .   From  these  times  t   and  t.  on  the 
3  SG     SR  a        b 

observed  wave  record,  one  moves  outward  from  the  region  of 
high  waves  to  the  first  zero-upcrossing  encountered.   This  in- 
sures an  integral  number  of  waves  in  the  group.   The  times  • 
corresponding  to  these  upcrossings  thus  mark  the  beginning 
and  end  of  the  group.   In  Figure  3,  the  group  begins  at  time 
t  ,  ends  at  time  t  ,  and  contains  five  waves. 

2 .   Selection  of  Identification  Window  Width,  r 

The  final  step  to  be  taken  in  applying  this  method  of 
wave  group  determination  is  the  selection  of  a  value  for  r, 
the  width  of  the  window  to  be  used  to  compute  the  short-term 
running  significant  height  (or  variance)  of  the  wave  record. 
It  was  considered  desirable  that  two  groups  containing  waves 
of  identical  heights  possess  identical  values  of  H 
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regardless  of  the  period  of  the  individual  waves.  The  window 
width  should  not  allow  a  wave  to  belong  to  more  than  one  wave 
group,  and  should  not  allow  a  single  large  wave  to  be  identi- 
fied as  a  group.  The  window  must  also  be  able  to  produce  a 
group  waveform  that  will  identify  two  distinctive  groups  that 
may  be  separated  by  only  one  wave. 

Seventeen-minute  wave  records  with  a  constant  wave 
period  were  artificially  produced  to  observe  the  effects  of 
differing  windows  on  constant  period  waves.   To  simulate 
groups,  the  amplitudes  of  selected  individual  waves  were  al- 
lowed to  vary.   It  was  determined  by  experimentation  with 
these  artificially  constructed  monochromatic  wave  records 

that  if  r  is  fixed  for  all  records,  H_„  differences  will  exist 

SG 

for  groups  differing  only  in  period.   This  is  to  say  that  two 

monochromatic  wave  groups  having  the  same  distribution  and 

sequence  of  wave  heights,  but  differing  only  in  their  periods, 

give  different  H    values  if  the  same  window  is  used.   How- 
3  SG 

ever,  identical  H  „  values  are  obtained  if  the  ratio  of  the 

SG 

width  of  the  window  to  the  period  of  the  waves  in  the  group 
remains  the  same.   It  was  therefore  concluded  that  the  window 
should  be  dependent  upon  the  period  of  the  waves  in  the  group. 

Normally,  the  average  period  of  the  waves  in  a  group 
occurring  in  a  17-20  minute  record  differs  from  one  group  to 
another.   Accordingly,  it  was  decided  to  choose  a  window 
equal  to  or  a  multiple  of  the  period  of  the  waves  in  the  domi- 
nant wave  groups,  which  represent  the  primary  concentration 
of  energy  in  the  record.   Since  past  studies  have  shown  that 
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the  average  period  of  the  waves  that  constitute  prominent 
groups  closely  corresponds  to  the  spectral  peak  period  of  the 
record  (Thompson,  1972;  Smith,  1974),  the  width  of  the  window 
(r)  was  determined  by  the  spectral  peak  period  of  the  wave 
record  T  . 

In  order  to  determine  whether  to  adopt  an  identifica- 
tion window  having  a  duration  equal  to  the  spectral  peak 
period  of  the  record,  or  some  fraction  or  multiple  thereof, 
a  series  of  experiments  using  different  window  lengths  was 
performed  using  the  same  artificial  records.   It  was  determined 
that  if  windows  of  different  duration  are  used  in  computing 
the  running  short-term  significant  height  in  a  given  record, 

different  H    values  result  for  the  same  group.   Observing 
S  G 

this  change  in  H   ,  and  since  one  can  expect  some  deviation 

S  G 

of  the  wave  periods  in  a  wave  record  from  the  spectral  peak 

period,  it  was  desired  to  choose  a  small  range  of  window 

widths  within  which  the  computed  values  of  H„   would  not 

SG 

greatly  alter  the  basic  shape  of  the  group  waveform.   The 
boundaries  of  the  individual  groups  should  not  be  affected. 

Several  trials  in  this  study  have  shown  that  a  small 
window,  one  that  approximates  one  wave  period  or  less,  may 
allow  the  same  wave  to  belong  to  two  groups,  which  was  con- 
sidered unacceptable.   It  was  also  found  that  with  a  small  win- 
dow, a  single  large  wave  may  be  identified  as  a  wave  group  in 
the  group  waveform.   An  advantage  with  a  shorter  identifica- 
tion window,  however,  is  that  the  group  waveform  is  more  re- 
sponsive to  period  changes  in  the  observed  wave  record.   It 
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was  found,  on  the  other  hand,  that  the  use  of  a  smoothing  win- 
dow which  approximates  three  wave  periods  or  greater  produces 
a  group  waveform  that  is  not  responsive  enough  to  period 
changes,  and  also  allows  the  possibility  that  two  wave  groups, 
which  are  distinctive  to  the  eye,  may  be  consolidated  into 
one.   A  wider  window  may  also  result  in  variance  values  that 
are  very  small,  and  this  can  eliminate  groups  of  relatively 
smaller  waves. 

For  this  study,  the  width  of  the  identification  window 
(r)  was  arbitrarily  chosen  to  approximate  twice  the  spectral 
peak  period  (a  compromise  between  extreme  windows)  to  the 
closest  second.   As  may  be  seen  in  Figure  4,  the  average  of 

the  differences  between  H  „  and  H  „  for  five  selected  wave 

SG        SR 

groups  decreases  nearly  linearly  for  windows  1  1/2  to  2  1/2 

times  the  peak  period.   Figure  4  presents  characteristics  of 

an  artificial  wave  record,  but  analysis  of  actual  records 

shows  similar  results.   Some  groups  show  linear  trends  of 

average  H   -H    for  more  extreme  windows;  however,  1  1/2  to 
S  G   S  R 

2  1/2  times  the  peak  spectral  period  appears  to  be  a  minimum 
range  of  r  for  linearity.   This  range  avoids  the  exponential 
characteristics  of  narrower  windows  and  the  flatter  non-re- 
sponsive wider  windows.   It  was  felt  that  windows  approximat- 
ing twice  the  peak  period  would  produce  a  satisfactory  group 
waveform  for  study. 

Two  additional  restrictions  were  placed  upon  the 
definition  of  a  group.   The  first  is  that  any  positive  identi- 
fication regions  (H    >  H   )  that  are  less  than  one-half  the 
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selected  identification  window  length  apart  will  be  treated 
as  one  group.   This  restriction  was  established  to  prevent 
the  possibility  of  one  wave  belonging  to  two  groups.   The 
final  restriction  is  the  requirement  that  a  wave  group  must 
contain  a  minimum  of  two  waves. 

To  summarize,  the  wave  groups  in  a  record  are  deter- 
mined by  the  following  steps: 

a.  Determine  the  mean  water  level  and  variance  of  the 
observed  waveform. 

b.  Determine  the  spectral  peak  period  and  select  an 
identification  window  to  approximate  twice  this  value, 

c.  Operate  the  identification  window  upon  the  observed 
waveform  to  produce  a  short-term  variance,  or  corre- 
sponding height  parameter.   This  value  (H    in  this 

S  G 

study)  is  plotted  at  the  midpoint  of  the  window, 

along  with  the  value  for  the  entire  record  (H    in 

S  R 

this  study).   The  identification  window  operates  from 
the  beginning  to  the  end  of  the  observed  waveform 
(the  record),  and  produces  the  group  waveform. 

d.  Regions  where  H    >  H    denote  possible  wave  groups. 

S  G      S  R 

Group  boundaries  are  defined  by  movement  outward  from 

the  times  of  intersection  of  H    and  H    until  a  zero- 

S  G        SR 

upcrossing  is  reached  in  the  observed  waveform. 

e.  Positive  identification  regions  must  be  separated  by 
more  than  one-half  the  identification  window  to  be 
considered  separate  groups,  otherwise  a  single  group 
is  presumed  to  be  present. 
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f.   Groups  must  contain  at  least  two  waves. 

Graphical  examples  of  the  final  product  of  this  method 
are  shown  in  Figures  5A-5E.for  records  with  different  spectral 
peak  periods.  All  of  the  above  steps,  including  restrictions, 
were  handled  by  computer  using  digitized  wave  data. 
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Ill  .   WAVE  RETCORD  AND  WAVE  GROUP  MEASURES 

Various  measures  of  the  characteristics  of  wave  records 
and  wave  groups  were  utilized  in  this  study,  and  these  are 
introduced  below.   Definitions  introduced  by  previous  authors 
or  in  common  use  are  used  whenever  possible. 

A.  WAVE    RECORD    MEASURES 

1 .  Spectral    Peak    Period,    T 

The  spectral  peak  period  is  the  period  that  repre- 
sents the  primary  concentration  of  energy  in  the  record.   The 

value  of  T   is  found  from  the  spectral  analysis  of  a  record. 

R  z 

2.  Variance,  V„ ,  and  Significant  Height,  H 

The  variance  of  a  wave  record  is  determined  with  the 

use  of  Equation  (1)  and  is  a  measure  of  the  wave  energy.  This 

value  is  used  in  the  determination  of  the  significant  height 

of    the    record,    H„„    =    4/v~~. 

SR       R 

3 .  Wave  Steepness,  y 

The  wave  steepness  measure  combines  the  effects  of 

T   and  V  ,  and  as  used  here  is  defined  by 
R         R 

27T  HSR 
TR 

B.  WAVE  GROUP   MEASURES 

1 .   Wave  Group  Duration,  D 

Wave  group  duration  is  the  time  interval  from  the 
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beginning  to  the  end  of  a  wave  group  measured  between  zero  up- 
crossings,  and  always  includes  a  whole  number  of  two  or  more 
waves.   This  is  illustrated  in  Figure  6. 

2 .  Number  of  Waves  in  a  Group,  N 

N  is  the  number  of  waves  in  a  group  of  duration  D, 
and  is  always  an  integer  of  two  or  more. 

3.  Wave  Group  Period,  T_ 
* Q 

T   is  the  average  period  of  the  individual  waves  that 
G 

constitute  a  wave  group:   T   =  D/N . 

G 

4 .  Height  of  the  Largest  Individual  Wave  in  a  Group, 

H 
max 


The  height  of  an  individual  wave  is  defined  as  the 

vertical  separation  between  a  crest  elevation  and  the  average 

elevation  of  the  two  adjacent  troughs.   H     is  the  maximum 

max 

height  occurring  in  a  group,  as  shown  in  Figure  6.   H     of  a 

3  '       3    *  max 

group  is  normalized  by  use  of  the  ratio  H    /HCn* 

max   S  R 

5 .   Average  Variance  Ratio,  AVR 

A  measure  to  represent  the  ratio  of  the  average  vari- 
ance of  a  wave  group  to  the  variance  of  the  record  as  a  whole 
was  needed.   Its  definition  and  determination  can  best  be  de- 
scribed with  the  aid  of  Figure  6  and  Equation  (4) .   A  dimen- 
sionless  number  is  formed  termed  the  average  variance  ratio, 

AVR.   This  ratio  may  be  found  by  squaring  and  averaging  the 

2 
H  „  values  from  t   to  t   and  relating  that  to  (H„ „)  ,  since 
SG  A      B  v  SR 

a  height  parameter  squared  differs  from  the  variance  by  only 
a  constant,  and  thus  is  proportional  to  the  energy.  Because 
of  ease  of  graphical  presentation  of  data,  a  measure  related 
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to  the  AVR,  termed  AVR  ,  will  be  used  as  shown  in  Equation  (5) 

% 

AVR   may  be  thought  of  as  the  percentage  by  which  the  average 

% 

variance  of  a  group,  V,  differs  from  the  variance  of  the 

G 

entire  record,  V  . 

a. 

This  procedure  may  be  expressed  as  follows: 


AVR  = 


V. 


VT 


tB/At 

£  2     (HSG)^ 

■ 

VAt 

,    where    k    = 

B             A 

2 
(HSR> 

At 

(4; 


V   -  V 
AVR„  =  -~; x  100  =  (AVR-1)  x  100 


V 
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IV.   WAVE  DATA  SELECTION 

The  wave  records  used  in  this  study  were  obtained  from 
the  Coastal  Engineering  Data  Network  (CEDN)  sponsored  by  the 
California  Department  of  Navigation  and  Ocean  Development  and 
the  University  of  California  Sea  Grant  Program.   Scripps  In- 
stitution of  Oceanography  provides  overall  direction  for  the 
program.   Four  network  stations  in  Southern  California  were 
utilized  and  their  locations  are  shown  in  Figure  7. 

Wave  records  were  obtained  with  bottom-mounted  pressure 
sensors  placed  at  depths  of  approximately  10  meters.   The 
sensors  do  not  possess  high-pass  filters  and  therefore  record 

0 

all  long  waves  present,  including  the  astronomical  tides. 
These  long-term  sea-level  variations  were  considered  to  be  a 
second-order  effect  when  compared  with  the  wave  heights  oc- 
curring and  were  neglected  in  this  study.   Wave  data  from 
these  sensors  are  digitized  by  CEDN  on  magnetic  tape  at  a 
sampling  rate  of  one  second,  and  are  calibrated  in  units  of 
pressure  (cm  of  water) .   Each  individual  record  was  of  approxi- 
mately 17-minutes  duration  (records  range  from  1024  to  1078 
seconds)  with  a  record  generally  taken  every  10  hours  at 
each  station. 

CEDN  produces  a  spectral  analysis  of  the  pressure  record, 
corrected  for  hydrodynamic  damping,  and  publishes  both  graphi- 
cal representations  (Figure  8)  and  tabular  values  (Figure  9) 
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of  the  energy  spectra  for  each  station  by  months.   The  wave 
records  used  for  wave  group  analysis  in  this  study,  selected 
from  these  CEDN  presentations,  were  chosen  for  the  display  of 
sharp  unimodel  spectra.   Individual  records  were  selected 
when  greater  than  30%  of  the  energy  occurred  in  one  CEDN 
period  band.   The  30%  limit  was  arbitrarily  chosen  for  the 
purpose  of  insuring  that  a  prominent  spectral  peak  existed 

within  the  record.  The  spectral  peak  period  used  in  the  deter- 
mination of  the  identification  window  width  was  chosen  by 
interpolation  of  the  CEDN  tabular  data  to  the  nearest  whole 
second. 

Computations  of  variance,  in  this  study,  were  derived 
from  the  pressure  records  provided  by  CEDN  and  were  not  modi- 
fied for  hydrodynamic  ^damping.   It  was  assumed  that  the 
dominant  period  determined  from  spectral  analysis  by  CEDN 
would  closely  correspond  to  the  dominant  period  of  the  pressure 
waveform  over  the  range  of  periods  dealt  with  (7  to  20  seconds) . 

Two  hundred  eight  records  were  selected  to  be  examined, 
which  resulted  in  the  identification  of  1643  individual  wave 
groups.   The  height-period  characteristics  of  these  records 
are  shown  in  Figure  10.   The  graph  also  contains  steepness 
curves,  where  wave  steepness  is  defined  by 


Y  =  *=  — 

l   g 


H. 


SR 


2  ' 


where   H    is  the  significant  height  and  T   is  the  spectral 


SR 


R 


peak  period.   The  steepness  of  rapidly  growing  and  fully 
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arisen  seas  in  the  range  of  periods  dealt  with  in  this  study 
is  greater  than  or  equal  to  0.02  according  to  wave  forecast- 
ing graphs  (Bretschneider ,  1958).   Accordingly,  the  wave 
records  examined,  when  corrected  for  hydrodynamic  damping, 
cover  a  wide  range  of  steepnesses,  from  fully  arisen  seas  to 
old  swell. 
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V.   ANALYSIS  AND  INTERPRETATION  OF  WAVE  DATA 

Each  record  was  examined  for  relationships  between  the 
group  and  record  measures  defined  in  Section  III,  and  for  in- 
terrelationships between  the  group  measures  themselves. 

Since  values  of  V   and  T   can  cover  a  wide  range,  these 

R         R 

record  measures  were  grouped  into  categories  of  variance  and 
period  for  statistical  handling.   The  records  were  first 
divided  into  five  variance  bands  as  shown  in  Figure  11.   Rec- 
ords with  H  „  greater  than  one  standard  deviation  from  the 
SR 

mean  significant  height  of  all  the  records  are  included  in 

2 
bands  where  the  variance  is  less  than  200  cm   or  greater  than 

2 
700  cm  .   The  remainder  of  the  records  were  divided  into  three 

2  2 

variance  bands:  250  to  350  cm  ,  400  to  500  cm  ,  and  550  to 

2  2 

650  cm  .   A  50  cm   separation  was  used  to  isolate  the  five 

bands  from  each  other.   These  five  intervals  are  equivalent 
to  significant  height  intervals  of  less  than  57  cm,  63  to  75 
cm,  80  to  89  cm,  94  to  102  cm,  and  greater  than  106  cm  respec- 
tively. 

Spectral  analysis  of  the  wave  records  performed  by  CEDN 

yields  1/T   to  an  accuracy  of  Af  =  +  0.0034  Hz;  thus,  over  the 

R  ~— 

range  of  wave  periods  of  interest,  7  to  20  seconds,  the  pre- 
cision with  which  the  periods  can  be  obtained  is  approximately 
+_  0.17  sec  and  +_  1.4  sec,  respectively.   It  would  be  desirable 

to  examine  the  relationship  of  discrete  values  of  T   with, 

R 
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or  incorporated  in,  other  parameters  (e.g.,  y  or  T _/T  ) ,  but 

G   R 

because  of  the  low  precision  with  which  T  .  is  known  (especial- 

R 

ly  at  longer  periods)  it  was  decided  to  deal  with  bands  of 

T  .   The  bands  were  chosen  to  follow  the  CEDN  produced  period 
R 

bands,  which  include  the  lower  limit  of  the  band,  but  not  the 
upper  limit.   The  distribution  of  records  with  respect  to 
these  period  bands  is  shown  in  Figure  12.   T   period  bands 
less  than  six  seconds  were  not  examined  because  of  the  effect 
of  extreme  hydrodynamic  damping.   There  was  found  to  be  insuf- 
ficient data  available  in  the  10-12  second  period  band  for 
statistical  treatment,  and  this  band  is  omitted.   The  6-8 
second  band  and  the  8-10  second  band  were  combined  because  of 
sparsity  of  data. 

A.   RELATIONSHIPS  OF  GROUP  MEASURES  TO  RECORD  MEASURES 
1 .   Average  Wave  Group  Period,  T 

Figure  13  shows  cumulative  distributions  (in  percent- 
age of  groups)  of  the  average  period  of  wave  groups  occurring 
in  records  having  spectral  peak  periods  falling  in  the  five 
indicated  bands.   It  appears  that  more  groups  have  periods 
falling  within  the  spectral  peak  period  band  than  in  any  of 
the  adjacent  bands.   However,  a  wide  spread  of  group  periods 
is  apparent,  such  that  the  longer  spectral  peak  period  bands 
do  not  possess  a  majority  of  the  total  groups.   The  spread  may 
be  caused  by  individual  groups  of  shorter  periodicity  produced 
toward  the  short-period  tail  of  the  frequency  spectrum,  and 
possibly  also  by  the  presence  of  secondary  wave  trains. 
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An  upper  cutoff  in  the  distribution  is  evident  for 
each  spectral  peak  period  band.   This  cutoff  closely  corre- 
sponds to  the  upper  limit  of  the  associated  peak  period 
band.   This  agrees  with  the  dispersion  concept  in  which  waves 
of  longer  periods  travel  with  greater  celerity,  and  thus 
groups  of  these  longer  periods  have  already  passed  the  wave 
gage.   The  distributions  shown  in  the  figure  appear  to  contra- 
dict observations  by  other  investigators  (Thompson,  1972; 
Smith,  1974)  that  the  average  group  period  is  equal  to  the 
period  of  maximum  energy  density.   An  explanation  may  lie  in 
the  selectivity  of  the  groups,  in  that  the  method  used  in  this 
study  may  identify  groups  that  were  not  significant  to  pre- 
vious investigators. 

The  relationship  between  T   and  the  variance  of  the 
wave  record  was  not  examined  because  a  wide  range  of  combina- 
tions of  height  and  period  are  possible  in  ocean  swell. 
2 .   Number  of  Waves  in  a  Group,  N 

The  number  of  waves  in  a  group  is  shown  in  the  cumula- 
tive percentage  curve  of  Figure  14  representing  the  208  wave 
records  examined.   The  figure  shows  that  the  modal  number  of 
waves  per  group  is  4  or  5,  and  that  10%  of  the  wave  groups 
have  9  or  more  waves.   The  maximum  number  of  waves  in  a  group 
was  found  to  be  26. 

Wave  records  with  similar  peak  period  bands  were 
separated,  and  the  cumulative  distribution  for  N  for  each 
band  shown  in  Figure  15.   The  distributions  of  N  for  each 
period  band  are  statistically  similar,  which  indicates  that 
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there  is  little  dependence  of  N  upon  the  peak  period.   There 
is  a  slight  tendency  for  longer  period  bands  to  have  a  greater 
percentage  of  wave  groups  of  larger  N. 

The  number  of  waves  per  group  as  a  function  of  the 
selected  variance  bands  is  shown  in  cumulative-distribution 
form  in  Figure  16.   As  may  be  seen,  there  appears  to  be  even 
less  variability  in  the  frequency  distribution  of  N  with  vari- 
ance than  there  is  with  spectral  peak  period.   It  is  concluded, 
therefore,  that  the  number  of  waves  in  a  group  is  essentially 
independent  of  both  the  dominant  period  and  the  energy  content 
of  a  wave  record. 

3 .   Amount  of  Time  Wave  Groups  are  Present,  ID 

The  cumulative  time  during  which  wave  groups  are  pres- 
ent in  a  record,  ED,  equals  the  combined  duration  of  the  wave 
groups ,  given  by 


:d  =  2      (N 


'G  i 


i=l 


where  n  is  the  number  of  wave  groups,  including  groups  trun- 
cated by  the  ends  of  the  record.   The  cumulative  distribution 
of  ZD  is  shown  in  Figure  17;  the  average  duration  of  the 
analyzed  wave  records  is  approximately  1075  seconds.   An  aver- 
age time  of  group  occurrence  was  found  to  approximate  one-half 
the  record  length,  with  a  maximum  deviation  of  about  +  200 
seconds . 

It  is  speculated  that  the  amount  of  time  groups  are 
present  does  not  depend  upon  the  variance  of  the  record,  since 
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N  and  T   are  not  dependent  upon  the  variance.   With  regard 
G 

to  the  relationship  of  ED  to  the  spectral  peak  period,  if  T 

were  to  increase,  longer  T   would  be  possible  and  this  would 

G 

tend  to  increase  the  duration  of  the  individual  groups.   How- 
ever, an  increase  in  T   would  tend  to  reduce  the  number  of 

R. 

waves  possible  in  a  17-minute  record,  and  therefore  the  number 

of  groups  possible.   It  thus  appears  that  a  relative  balance 

occurs  between  the  individual  group  duration  and  the  number 

of  groups  per  record  such  that  a  change  in  T   would  be  expected 

R 

to  have  little  effect  upon  Zd. 

4 .   Average  Variance  Ratio  Percentage,  AVR 

As  stated  above,  AVR   is  a  measure  of  the  energy  in 

% 

an  individual  wave  group  relative  to  the  energy  in  the  record; 

positive/negative  values  indicate  that  the  wave  group  has  more/ 

less  energy  than  the  wave  record.   The  distribution  of  AVR 

for  all  1643  groups  identified  is  shown  in  Figure  18.   It  is 

seen  that  approximately  84%  of  the  groups  have  an  energy  level 

greater  than  the  record  as  a  whole. 

The  groups  were  once  again  separated  into  peak  period 

bands,  and  the  distribution  of  AVR   by  these  bands  is  shown 

% 

in  Figure  19.   It  may  be  seen  that  the  AVR   distributions  are 

% 

similar,  indicating  that  AVR   is  independent  of  the  spectral 
peak  period  of  the  record. 

Figure  20,  on  the  other  hand,  shows  a  distinctive  rela- 
tionship between  the  AVRo  occurrence  and  the  variance  of  the 

wave  record.   It  may  be  seen  that  the  greater  the  value  of  V  , 

R 

the  more  heavily  weighted  the  distribution  of  groups  is  toward 


32 


higher  AVR   values.   Stated  in  more  practical  terms,  the 
greater  energy  a  wave  train  possesses,  the  greater  is  the 
possibility  that  individual  wave  groups  will  have  propor- 
tionately higher  energy. 

5.   Maximum  Wave  Height  in  a  Group,  H 

2 C£ —  max 

The  value  of  H    ,  the  single  highest  wave  in  a  group, 
max  F 

may  be  expected  to  be  statistically  dependent  upon  the  vari- 
ance of  the  record  from  consideration  of  the  Rayleigh  dis- 
tribution applied  to  ocean  waves.   Thus,  as  variance  increases 
the  probable  maximum  height  obtainable  in  the  record  should 
be  expected  to  increase  with  the  square  root  of  the  variance. 
This  relationship  is  generally  indicated  in  Figure  21  for  the 
upper  and  lower  limits  of  the  distribution.   A  change  in  the 
spectral  peak  period  for  the  same  number  of  waves  would  not 

be  expected  to  affect  the  H     values  if  the  energy  level  re- 

max 

mained  the  same  in  the  record.   This  is  inferred  from  theoreti- 
cal considerations  and  is  not  presented  graphically. 

B.   INTERRELATIONSHIPS  OF  GROUP  MEASURES 

An  examination  of  the  six  possible  interrelationships  of 

the  four  group   measures  AVR„ ,  T  ,  N,  and  H     was  conducted. 

3    e  %    G  max 

The  principal  observations  follow. 
1.   AVR„  and  T 


Figures  22-26  depict  the  relationship  of  AVR   to  the 

% 

average  group  period  for  the  spectral  peak  period  bands  indi- 
cated.  It  appears  that  the  higher  values  of  AVRo  generally 

coincide  with  the  groups  whose  T   approaches  the  spectral  peak 

G 

period.   The  figures  also  show  a  relatively  sharp  period 
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cutoff  at  the  upper  end  of  the  spectral  peak  period  band, 
above  which  few  groups  and  low  AVR   values  occur. 

2.  AVR„  and  N 

Figures  27-33  show  the  relative  distribution  of  AVR 
as  a  function  of  the  number  of  groups  for  a  constant  N.   It 
may  be  seen  that  as  the  number  of  waves  in  a  group  increases, 
the  distribution  migrates  toward  higher  AVR   values  and  be- 
comes  less  peaked. 

It  has  been  shown  that  AVR0  is  highly  dependent  upon 

the  variance  (Figure  20) ;  for  an  increase  in  variance  there 

must  be  a  corresponding  increase  in  the  percentage  of  groups 

with  high  AVR  .   This  redistribution  of  groups  with  respect 
% 

to  AVR   must  occur  with  no  change  in  the  distribution  with 
% 

respect  to  N,  since  it  has  also  been  shown  that  N  appears  to 
be  independent  of  the  variance  (Figure  16).   This  is  evident- 
ly accomplished  by  the  shifting  of  groups  having  a  particular 
N  to  higher  AVR   values  within  the  ranges  of  AVR0  implied  by 
Figures  27-33.   This  is  illustrated  by  comparison  of  the  dis- 
tributions shown  in  Figures  34  and  35  for  the  extreme  vari- 

2  2 

ance  bands  (V   <  200  cm   and  V   >  700  cm  ) .   The  figures 
R  R 

give  the  approximate  percentage  of  groups  in  a  number  of  AVR 

bands,  and  show  that  for  an  increase  in  the  total  energy  of  a 
record,  a  much  wider  range  of  relative  wave  group  energy 

(AVRo )  is  possible.   It  may  also  be  noted  that  there  are  ap- 

-3 

parent  limits  to  AVR   for  each  individual  value  of  N. 

% 

3.  AVRa  and  H    /E„„ 

% max   SR 

This  relationship  is  shown  in  Figure  36.   Separated 
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bands  of  AVR   (-0.25  to  0,  0.50  to  0.75,  1.00  to  1.25,  2.00 

to  2.25,  and  >  3.00)  were  chosen  to  emphasize  the  vertical 

spread  and  trend  in  H    /Hnn  as  AVR   increases.   The  figure 
e  max   SR        % 

shows  that  groups  with  large  AVR   are  more  likely  to  have  an 

extreme  H     than  groups  with  small  AVR„  .   In  addition,  a 
max  % 

particular  value  of  H    /H„_  cannot  be  used  as  a  measure  of 
e  max   SR 

the  relative  energy  in  a  group,  but  may  suggest  a  minimum 

energy  level  possible.   It  may  also  be  noted  that  the  maximum 

height  in  a  group  may  be  less  than  the  significant  height  of 

the  record;  the  figure  shows  that  H     is  dominantly  less 

max 

than  H    for  AVR   values  less  than  1.00,  but  greater  than 
S  R  % 

H  „  for  AVR„  values  of  2.00  or  greater. 
SR         % 

4 .  N  and  T 

~~ — ^— — — ~~  Qj 

Figures  37-41  show  a  general  trend  that  the  number  of 
waves  in  a  group  will  probably  be  greater  when  the  average 
period  of  the  group  approaches  the  spectral  peak  period  of 
the  record.   This  is  a  comfortable  conclusion  since  the  far- 
ther away  the  average  period  of  a  group  is  from  the  spectral 
peak  period,  the  fewer  are   the  number  of  waves  that  would  be 
anticipated  from  wave  interference  considerations.   Wave 
groups  with  average  group  periods  greater  than  the  spectral 
peak  period  band  are  once  again  seen  to  be  few  in  number  and 
have  small  N. 

5.  w  and  H    /H__ 

max   SR 

The  more  waves  in  an  individual  wave  group,  the  greater 
is  the  chance  of  obtaining  an  extreme  wave  height.   This  is 
shown  in  Figure  42.   For  a  given  N,  there  is  an  apparent  lower 
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limit  to  the  value  of  H    /H__,  and  for  a  given  H    /H__  a 

max   SR  max   SR 

wide  range  of  N  is  possible.   It  may  also  be  noted  that  for 

groups  with  12  or  more  waves,  the  maximum  height  in  the  group 

is  likely  to  be  greater  than  the  significant  height  of  the 

record. 

6.   T   and  H    /HOT, 
G max   SR 

The  final  interrelationship  was  examined  without  the 

aid  of  a  figure.   H   • /H„„  was  shown  to  increase  with  large 
3        max   SR  * 

AVR   (Figure  36)  and  large  N  (Figure  42).   However,  AVR   and 
%  % 

N  were  shown  to  be  independent  of  the  spectral  peak  period 
(Figures  19  and  15),  but  dependent  upon  how  closely  the  aver- 
age group  period  approximates  the  spectral  peak  period 
(Figures  22-26  and  Figures  37-41).   Therefore,  extreme  rela- 
tive wave  heights  (i.e.,  large  values  of  H    /H„„)  should  be 
r  max   SR 

expected  when  the  group  period  approaches  the  spectral  peak 
period,  regardless  of  the  value  the  spectral  peak  period  as- 
sumes . 

C.   OTHER  RELATIONSHIPS 

In  order  to  determine  how  wave-group  properties  vary  over 
the  life  of  an  arriving  swell  train,  a  time-series  examina- 
tion was  made  of  12  synoptic  swell  trains,  each  sampled  by  a 
17-rainute  record  at  approximately  10-hour  intervals.   Each 
train  was  selected  on  the  basis  of  displaying  a  well-defined 

wave  height  maxima  and  a  systematic   decrease  of  T   with  time. 

R 

It  was  found  that  as  the  wave-record  measures  of  a  synoptic 
wave  train  change,  the  associated  wave-group  properties  change 
in  a  manner  that  would  be  expected  from  consideration  of  the 
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analysis  described  above.   The  presentation  of  this  work 
is  not  included  since  no  additional  information  concerning 
wave-group  properties  was  revealed. 

A  comparison  of  records  with  respect  to  the  four  differ- 
ent wave  gage  sites  likewise  resulted  in  the  conclusion  that 
wave-group  relationships  are  independent  of  location.   Differ- 
ences in  degree  and  direction  of  site  exposure  to  open  ocean 
wave  conditions  resulted  mainly  in  wave  steepness  differences, 
which  influence  wave  groups  only  as  the  height  parameter 
changes . 
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VI.   SUMMARY 

A  method  of  digitally  analyzing  an  ocean  wave  record  for 
wave-group  characteristics  with  the  use  of  a  moving  group 
identification  window  was  developed.   The  window  was  chosen 
to  be  dependent  upon  the  spectral  peak  period  of  the  record. 
The  short-term  variance  of  the  record  within  this  moving  win- 
dow was  compared  to  the  variance  of  the  entire  record  to 
identify  intervals  of  apparent  energy  concentration.   These 
intervals  were  considered  to  be  related  to  the  occurrence  of 
wave  groups . 

Two  hundred  eight  wave  records  that  exhibited  unimodel 
spectra  were  selected  for  analysis.   The  recording  sites  were 
four  stations  on  the  Southern  California  coast  that  are  part 
of  the  Coastal  Engineering  Data  Network  system  of  the  State 
of  California.   Wave  record  measures  (the  spectral  peak  period 
and  the  variance),  as  well  as  wave  group  measures  (the  average 
group  period,  the  number  of  waves,  the  maximum  wave  height, 
and  the  relative  group  energy)  were  defined,  determined  by 
computer  analysis,  and  then  related  to  one  another. 

It  was  found  that  the  number  of  waves  in  a  group  is  inde- 
pendent of  both  the  spectral  peak  period  and  the  variance  of 
the  wave  record.   However,  the  amount  of  energy  contained  in 
wave  groups  relative  to  that  in  the  record  increases  as  the 
total  energy  of  the  record  increases.   Also,  average  group 
periods  of"  greater  than  the  spectral  peak  period  are  unlikely 
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to  occur,  but  group  periods  less  than  the  peak  period  are 
not  uncommon.   Both  an  increase  in  group  energy  relative  to 
the  wave  record  and  also  in  the  number  of  waves  in  a  group 
increase  the  possibility  of  obtaining  an  extreme  wave  height 
in  a  group  relative  to  the  significant  height  of  the  record. 
Further,  the  closer  that  the  average  wave  group  period 
approximates  the  spectral  peak  period  of  the  record,  the 
higher  is  the  relative  energy  and  number  of  waves  possible 
in  the  group.   Further  study  of  this  relationship  is  extreme- 
ly desirable  since  it  is  evident  that  all  wave  group 
characteristics  are  dependent  upon  it. 

In  view  of  these  findings,  it  would  therefore  be  expected 
that  in  a  synoptic  wave  train  the  largest  ratio  of  group  to 
record  energy,  the  largest  ratio  of  maximum  wave  height  within 
a  group  to  significant  height  of  the  record,  and.  the  largest 
number  of  waves  in  a  group  would  occur  at  the  time  of  peak 
energy  arrival.   This  conclusion  may  be  expected  to  apply 
regardless  of  the  spectral  peak  period. 
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Figure  2.   Group  Waveform  in  Relation  to  Observed  Waveform. 
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Figure  3.   Determination  of  Wave  Group  Boundaries 
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Figure  4.   Average  Differences  in  H    and  H    as  a  function 
of  Window  Width  (5  selected  wave  groups) . 
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Figure  5A.   Example  of  Group  Waveform  Analysis 
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Figure  5B.   Example  of  Group  Waveform  Analysis 
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Figure  5C.   Example  of  Group  Waveform  Analysis. 
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Figure  5D.   Example  of  Group  Waveform  Analysis. 
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Figure    5E.       Example    of    Group    Waveform    Analysis. 
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Figure  6.   Determination  of  AVR  and  H 


max 


51 


g 


< 


Figure  7 


Location  of  CEDN  Stations  (from  Seymour,  et  al , 
19771  . 
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WAVE  ENERGY  SPECTRA  DURING  AUGUST  1375 
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Figure  8.   Graphical  Presentation  of  CEDN  Spectra 
(from  Seymour,  et  al  ,  1977). 
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Tabular  Presentation  of  CEDN  Spectra,  Ocean  Beach, 
August  1-9,  1976  Cfrom  Seymour,  et  al ,  1977). 
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Figure  10.  Height-Period  Distribution  of  Selected  Wave  Records. 
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Figure  11.  Distribution  of  H    in  Wave  Records  Analyzed  and 
Variance  Bands  Used. 
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Figure    12. 


Distribution  of  Spectral  Peak  Period  in  Wave 
Records  Analyzed. 
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Figure  13.   Distribution  of  T   by  Spectral  Peak.  Period  Bands. 
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Figure    14.       Distribution    of    N     (1643    groups) 
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Figure  15.  Distribution  of  N  by  Spectral  Peak.  Period  Bands 
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Figure  16.   Distribution  of  N  by  Variance  Bands 
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Figure  17 
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Figure  18.   Distribution  of  AVR   (1643  groups) 
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Figure  19.   Distribution  of  AVR   by  Spectral  Peak.  Period 
Bands . 
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Figure  20.   Distribution  of  AVR   by  Variance 
Bands . 
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Figure    21.       Variance    va    H 
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Figure    22.       AVR„    vs    T_     (.18-22    sec    band) 
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Figure    23.       AVR      vs    T       C16-18    sec    band) 
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Figure  24.   AVR   vs  T   (.14-16  sec  band) 
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Figure    25.        AVR.     vs    T        C12-14    sec    band) 
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Figure    26.        AVR„     vs    T        (6-10    sec    band) 
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Figure  27.   Distribution  of  AVR   for  N  =  2 


72 


[ 


I 


o 


00 '09    00'Sh    00 'OS    00 'SI 

SdnoyQ  jo  H39wnN 


o 
o 

M     c 


o 
o 

D    • 

in 


o 
o 

* 


<\J> 


a 

o 


o 

o 


a 
a 

« 


00  *o' 


Figure    28.       Distribution    of    AVR      for    N    =    3 
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Figure  29.   Distribution  of  AVR   for  N  =  4. 
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Figure  30.   Distribution  of  AVR   for  N  =  5 
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Figure  31.   Distribution  of  AVR   for  N  =  6 
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Figure  32.   Distribution  of  AVR   for  N  =  8 
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Figure  33.   Distribution  of  AVR   for  N  =  10 
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Figure    34.       AVR      vs    N    for    V       <    200    cm    Cpercentage    x    10       ) 


79 


U"> 


in 


W 


in 


CD 


I/) 


CD 


CO 


in 


in 


m 


CM 


in 


in 


en 
cn 


in 


CM 
CM 

CO 

CM 
CM 

in 

CM> 

CO 

t-1 

CM 
CM 

CO 

CE 

CO 

cn 

CO 

O 

o 

m     • 

1-t 

CO 

1"* 

CM 
CM 

cn 

CO 

CM 
CM 

3" 

in 

cn 
cn 

o 

o 

—     • 

o 

CD 

3* 
3" 

00 

cn 

V* 

CM 

o 

o 

m 


(N)   dnoyo  y3d  S3aum 


CM 


CO 


in 


o 
o 

—     • 

CO 


o 
o 

in 


o 
o 


o 
o 

°cn 


o 

o 

LCM 


Figure  35.   AVR   vs  N  for  V   >  700  cm  (percentage  x    10   ) 
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Figure    36.       AVR      vs    H         /H 
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Figure  37.   N  vs  T   (18-22  sec  band) 
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Figure    38.       N    vs    T       C16-18    sec    band). 
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Figure    40.       N    vs    T       C12t-14    sec    bandl. 
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Figure    41.       N    vs    T       C6-1Q    sec    bandl . 
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